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Abstract 

Tailoring the properties of correlated oxides is accomplished by chemical doping, pressure, tem¬ 
perature or magnetic held. Photoexcitation is a valid alternative to reach out-of-equilibrium states 
otherwise inaccessible. Here, we quantitatively estimate the coupling between a lattice distortion 
and the charge-transfer excitation in (La 2 Cu 04 _|_ 5 ). We photoinduce a coherent La ion vibration 
and monitor the response of the optical constants in a broad energy range, providing quantitative 
information on the electron-phonon matrix element that can be compared to theoretical models. 
We propose the same methodology to probe electron-electron interactions in other materials. 
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The non-trivial interplay between low- and high-energy degrees of freedom is one of the 
most distinctive characteristics of strongly correlated systems HHII. An elegant way to inves¬ 
tigate this phenomenology is provided by pump-probe experiments. In Coherent Fluctuation 
Spectroscopy (CFS), a pump pulse induces a coherent excitation whose effect on the optical 
properties is monitored in a broad energy window via femtosecond (fs) reflectivity. Direct 
information on the coupling between the coherent mode and the electronic excitations in the 
accessible spectral range is obtained. Coherent modes that can be studied include charge 
fluctuations in coherent charge fluctuation spectroscopy (CCFS) and lattice vibrations in 
coherent lattice fluctuation spectroscopy (CLFS). For example, Mansart et al. [5] argue that, 
in order to obtain information on possible pairing actors, rather than using the traditional 
approach of changing a bosonic excitation and monitor the effect on the superconducting 
properties (as in the classical isotope effect [HHTnjh one can make use of the reciprocity prin¬ 
ciple and study how a perturbed condensate affects bosonic excitations. Thus, using CCFS, 
it was found that superconducting quasiparticles are coupled to a high-energy excitation on 
the Mott scale (2.6 eV). Although these experiments provided evidence for the coupling of 
the condensate quasiparticles to other excitations, the magnitude of the coupling has not 
been evaluated so far. 

Since the most common mechanism to create coherent excitations involves Raman matrix 
elements [nHi6], one can apply CFS to virtually any excitation which can be triggered by 
impulsive or displacive stimulated Raman scattering, i.e. any excitation that is Raman active 
and has a period longer than the available pump pulse duration. In impulsive stimulated 
Raman scattering (ISRS) [HI [12], the pump pulse can be modeled as a time-dependent 
electric held acting on the solid, E(t) = Re[£(t)e*®'^*/% where El/K corresponds to the 
central laser frequency of the pump and £{t) describes the shape of the 45-fs pump pulse 
used in our experiments. In the following, we use E to denote energies in the eV range and 
uj to indicate frequencies of coherent oscillations in the THz regime. In ISRS, the Raman 
interaction converts the incoming optical pulse into an impulsive force E{t) acting on the 
ions. Assuming, for simplicity, a single ion moving in the volume Uj, the force in transparent 
media averaged over one period of the excitation is PH IQ] 

F(t) = ■ £•(«). (1) 

where 2 ; is the ionic coordinate, e is the dielectric tensor of the material and its derivative is 
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FIG. 1. (a) Real part and (b) imaginary part of the dielectric function e from spectroscopic 

ellipsometry measurements. The vertical line indicates the pump energy, the shaded area denotes 
the probe spectrum. The dashed blue lines are fits to the data with the model of Ref. [20] . 


taken at the pump-laser energy El. The quantity |e ■ ■ e*p (e being a unit polarization 

vector) determines the cross section for spontaneous Raman scattering [TTl HH] and its abso¬ 
lute magnitude has been determined only in few cases [19] , with very demanding continuous 
wave experiments performed at many different laser frequencies. Here, we show that this 
information can be obtained from pump-probe spectroscopy in a single experiment. 

Assuming a single oscillator of mass M and frequency cuq, the impulsive force from Eq. Q 
produces an oscillation with maximum amplitude 
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where the fluence E measures the energy per unit surface arriving at the sample with each 
pulse. The dielectric tensor = e - e ■ e* in Eqs. Q and (|^ is projected according to the 
polarization of the pump. The ion movement produces oscillations of the reflectivity given 
by: 
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where Spr is the dielectric tensor projected according to the polarization of the probe and 
t is the time delay between pump and probe. The hrst factor on the right hand side of 
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Eq. (|^ is known from Fresnel’s equations. Taking polarizations such that £pr = £pu and 
evaluating Eq. ([^ at the time of the maximum amplitude and at the energy of the pump one 
can eliminate z from Eqs. (|^ and ([^ and obtain an equation for the squared magnitude of 
in terms of experimental quantities. Eq. ([^ then gives access to the full resonant 
Raman prohle. In practice, since relating the reflectivity to £ requires a Kramers-Kronig 
(KK) transform, it is more convenient to use a model dielectric function as we show below. 

To test these ideas, we performed ISRS experiments on a cuprate parent compound and 
we report a strategy to obtain the absolute coupling between low- and high-energy degrees 
of freedom using CLFS. An Aig La-phonon is coherently excited in the (slightly doped) 
cuprate parent compound (La 2 Cu 04 + 5 , or LCO) by using fs pump pulses polarized parallel 
or perpendicular to the Cu-0 planes. The consequent real-time oscillation of the optical 
reflectivity is monitored in the charge-transfer (CT) region to retrieve the energy and po¬ 
larization dependence of its strength. In this situation, the Raman matrix elements carry 
microscopic information on the electron-phonon matrix elements [12]. Currently available 
methods for evaluating the electron-phonon coupling allow averaged information over all 
lattice vibrations to be obtained [2IH23]. In contrast, we show that CLFS provides infor¬ 
mation on the electron-phonon coupling of specific phonons in a correlated system, giving 
strong anchor points to theory and providing a useful knowledge on the selectivity of the 
pump energy and polarization in the excitation of a mode. 

A slightly doped single crystal of LCO was grown as described in the SI, oriented in a 
Lane diffractometer, cut along a plane containing the a- and c-axis and polished to optical 
quality. The Neel temperature was determined to be Tjv = 260 K, which corresponds to a 
doping hole content of 6 ■ 10“^ or an oxygen content of 5 = 3 ■ 10“^ according to Ref. [21]. 
Static optical data between 0.5 eV and 4 eV with a resolution of 40 meV were obtained 
at room temperature using spectroscopic ellipsometry in a hxed polarizer geometry for 70° 
and 75° angles of incidence. Reflectivity measurements with 1.55 eV pump and white-light 
(1.7-2.8 eV) probe were performed at a temperature of 10 K and 6 KHz repetition-rate in a 
closed cycle liquid helium cryostat. The overall resolution of the experiment was 45 fs (see 
SI). 

Figure [T] shows the real and imaginary parts of the dielectric function e of the sample under 
study as measured by spectroscopic ellipsometry, revealing the strong optical anisotropy of 
the material. The in-plane imaginary part £2 of the dielectric function is dominated by 
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FIG. 2. (a) In-plane reflectivity change of La 2 Cu 04 _|_ 5 , shown in color-coding, (b) Time-dependent 
reflectivity change for selected probe energies. Traces are integrated over 0.15 eV wide slices 
centered on the horizontal lines in (a). The pump light is polarized parallel to the c-axis. (c) 
Oscillatory contribution of our data for E = 2.1 eV. (d) Eigenvector of the La Aig phonon 

mode identified as the oscillation in our data, after |25) . 


the oxygen 2p to a copper 3 dx' 2 -y 2 CT around 2.15 eV. Although not strictly of excitonic 
nature, Falck et al. m have described the CT excitation as creating strongly interacting 
electrons and holes because of the reduced width of the CT peak. The hnite width of 
the peak as well as a shoulder around 1.6 eV reveal the presence of excess oxygen in our 
sample, the actual composition being La 2 Cu 04 + 5 . Comparing 82 to static data by Uchida 
et al. [2S] we determine that our sample has an effective doping of less than a percent, 
consistent with the Neel temperature estimate. In contrast to the in-plane behavior, £2 
is rather featureless along the c-axis. As expected, the sample is not a good conductor 
in this direction, although it shows a feature at the energy of our pump beam. To our 
knowledge, there are no specihc transitions assigned to the c-axis spectrum at our pump 
energy, although the in-plane shoulder present at the same energy can be attributed to a d-d 
transition producing trapped carriers ^7]. Nevertheless, it is a general feature of cuprate 
superconductors that excitation along the c-axis is an effective way of exciting phonons due 
to enhanced electron-phonon coupling in this conhguration |28j . 

For the time-resolved study the sample was excited at 1.55 eV, indicated by the vertical 
line in Fig. The shaded area highlights the probe range of the experiment. We study 
the case of the pump polarized both in-plane and perpendicular to the plane. In the latter 
case incoherent excitations are suppressed due to the insulating character allowing a better 
observation of coherent oscillations, so we start by analyzing this case. Figure [^a) shows 
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the in-plane reflectivity changes of the sample at 10 K for = 7.8 mJ cm“^ incident fluence 
(polarized out-of-plane) as a function of time delay and probe energy for the first 2.5 ps 
after laser excitation. 

The dominant feature is a strong negative peak around 2.15 eV probe energy that decays 
within ~ 0.5 ps, very similar to the dynamics observed by Novelli et al. [1] in a comparable 
experiment. The energy maximum coincides with the maximum of 82 , but the line shape 
is much narrower, indicating that the CT feature in Fig. [^b) is indeed superimposed by 
the shoulder-like feature around 1.6 eV unrelated to the CT. The incoherent part of the 
response can be understood as being caused by the blocking of excitations across the CT 
peak due to the depletion of the initial state (population of the final state) by excitations 
along the c-axis to (from) the intragap states observed in the static data. The incoherent 
peak stretches into a short-lived tail towards higher probe energies (figure top), whereas 
towards lower energies two vertical lines indicate the presence of a coherent excitation. This 
excitation, which we identify as a coherent La Aig phonon (see SI) and which will be our 
primary interest in this work, is present over the whole probe spectrum, as can be seen from 
the time traces for fixed probe energy in Fig. [^b). 

We used a singular value decomposition (SVD) algorithm to separate the differential re¬ 
flectivity matrix into an incoherent part and an oscillatory part (Ai?/i?)°®'^ {t, E) as detailed 
in the SI. A damped oscillation was fitted to the dominant time traces in the SVD of the 
differential reflectivity matrix and the weight of each SVD trace was used to reconstruct the 
oscillatory part of the reflectivity. The reconstructed damped oscillation at a fixed probe 
energy is shown in Fig. |^c). The decay time is r = 170 fs and the natural frequency is 
(Uq = 4.8 THz, which agrees with the frequency reported in spontaneous Raman scattering 
for the Aig phonon in which the La atoms move along the c-axis (see SI), as depicted in 
Fig.gd). 

The energy profile of the oscillation extracted from the SVD is shown in Fig. Ha). It 
shows a resonance around 2.1 eV, which matches the CT feature seen in our static data. 
The energy profile for the same oscillation observed with the pump polarized in-plane is also 
shown. Due to the lower oscillation-to-background ratio, it can only be obtained at a few 
energy values, integrating over a wide energy range for each point. According to Eqs. ([^ 
and ([^, the absolute Raman profile can be obtained when pump and probe have the same 
polarization, but the profile shape Eq. (|^ depends only on the probe direction. Indeed, we 
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obtain very similar line shapes. Therefore, we rescale the prohle for the pump beam polarized 
along the c-axis to match the intensity of the case with the pump polarized in-plane, which 
allows the full differential reflectivity prohle to be obtained as if the pump was polarized in 
the plane. The next step is to hnd the Raman prohle ^e{E) = ^{E)zm- While in principle 
this can be extracted from the diherential rehectivity prohle using KK transformations, we 
hnd it more convenient to extract the information by htting the data to a simple model. As 
a bonus, this allows the Raman prohle to be related to the electron-phonon matrix elements. 

We model the in-plane static dielectric function e using an interacting polaron model as 
described by Falck et al. EDI- The absorption edge is modeled in terms of transitions between 
a valence and a conduction band modihed by Coulomb interactions which for simplicity are 
described by a momentum-independent matrix element yielding 


E^Uj) — Eoo + S ■ 


Gju;) 

l-UehG{uj)' 


(4) 


Here, the oscillator strength is parametrized by s = ine'^xlf^/vcu, where Xeh is a dipole 
matrix element between the Wannier orbitals of the two bands and Vcu is the volume per 
Cu atom. We also dehne the electron-hole Green’s function, 

Eg + W 

Eg 


where, due to the two-dimensionality of the system, we assumed a hat joint density of states 
D = 2/W conhned to a band of width W above the gap energy Eg. The bandwidth W of the 
absorption band is the sum of the electronic (single-particle) hole and electron bandwidths. 
The linewidth of the transition is given by half the (two-particle) phenomenological damping 
7 , which we assume has a simple linear dependence on the frequency, 7 = a -|- 6 ■ cu. 

Notice that rigorously speaking the model is for a stoichiometric compound, while we 
have a small oxygen surplus. However, the doping is so small that its ehect can be phe¬ 
nomenologically absorbed in small changes of the damping and other parameters. According 
to the Born-Oppenheimer approximation it is assumed that the electronic parameters of the 
model, like the density of states D and gap Eg, depend parametrically on the coherent ionic 
motion excited by the pump. 

We hnd that we can simultaneously ht the equilibrium optical response in Fig. and 
the differential rehectivity prohle at the maximum of the oscillation at time Im in Fig. |^a) 
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with only two parameter changes, 6Eg = 1.2 meV, 6D = 8.7 ■ 10“^ eV“^. Physically, these 
are also the parameters expected to be most sensitive to the La motion. Fig. |^b) shows 
the differential dielectric function A^i corresponding to the oscillation amplitude at Im- 
Neglecting the absorption part, we can represent its value at the pump energy as Ae^El) = 
^{El)zm = 4.4 X 10“^. This equation together with Eq. (|^ (using for M the mass of one 
La atom moving in the volume Vi = nc/4, Vc being the orthorhombic unit cell volume) yields 
a maximum amplitude zm = 3.2 x 10“^ A. We can use this value to calculate the absolute 
Raman prohle || = Ae{E)/ zm-, shown in the right scale in Fig. j^b). We are also able to 
evaluate the electron-phonon matrix elements dEg/dz = bEgjzu = 2.8 eVdescribing 
the change of the CT gap with the La motion and dD/dz = 6D/zm = 2.0 eV“^A“^ 
describing the change in density of states. 

To check if this estimate is reasonable we use a simple ionic model and assume dEg/dz 
is approximately given by the change of the difference in Madelung energy among planar 
Cu and O atoms when the La is moved in the x direction (see SI). We obtain dEg/dz = 
4.3 eVA-i , in reasonable agreement with the experimental value. We attribute the discrep¬ 
ancy to the strong covalency of the Cu-0 bond which will make the CT energy less sensitive 
to the change of Madelung energy with respect to what the simplihed ionic estimate suggests. 

To summarize, we have shown that coherent fluctuation spectroscopy can be used to 
determine absolute Raman prohles much easier than static Raman scattering. One potential 
application of the Raman prohle is the manipulation of the electronic ground state by direct 
excitation of specihc phonons [294133] . In the case of Raman excitation, the Raman prohle 
allows to hnd the pump laser energy at which the applied force is maximum, thus allowing to 
reduce huence and unwanted heating. For example. Fig. [^b) together with Eq. ([^ suggests 
that exciting at El = 2 eV the oscillations can be enhanced by an order of magnitude using 
the same huence E. 

We also show that the Raman prohle can be described via changes in few key-electronic 
parameters describing the equilibrium optical properties, CT energy and density of states. 
This allows the microscopic electron-phonon matrix elements to be estimated in a correlated 
system, hnding agreement with the upper limit given by a simple ionic model. As a bonus, 
we conhrm that ISRS is the mechanism at work on generating the ionic motions. We are 
not aware of previous works where it was checked for a strongly correlated system that the 




FIG. 3. (a) Differential reflectivity energy profile of the coherent oscillation at t = tM for out- 
of-plane and in-plane excitation, and fit. The data with the pump || c was rescaled to match the 
intensity of the data with the pump _L c. (b) Differential dielectric function Aei corresponding to 
in-plane pump and probe polarizations at t = tM (left scale). The absolute right scale was obtained 
as explained in the text. 

amplitude of the oscillations is correctly given by this mechanism as onr analysis implies. 
This is important given that other mechanisms have been proposed in the past [12]. Finally, 
onr resnlts pave the way to a systematic investigation of the interplay between charge |S| 
and/or structural fluctuations and electronic excitations in complex oxides. 

We acknowledge professor D. van der Marel for allowing us to use his optical spectroscopy 
equipment for the static experiments. Work at LUMES was supported by ERG starting grant 
USED258697 (F.C.) and the NCCR MUST, a research instrument of the Swiss National 
Science Foundation (SNSF). 
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SUPPORTING INFORMATION 


Sample preparation 

Polycrystalline La 2 Cu 04 samples were prepared by a solid state reaction. The starting 
materials La 203 and CuO with 99.99% purity were mixed and ground followed by a heat 
treatment in air at 900-1050°C for at least 70 hours with several intermediate grindings. 
The phase purity of the resulting compound was checked with a conventional x-ray diffrac¬ 
tometer. The resulting powder was hydrostatically pressed into rods (7 mm in diameter) 
and subsequently sintered at 1150°C for 20 hours. The crystal growth was carried out using 
an optical floating zone furnace (FZ-T-IOOOO-H-IV-VP-PC, Crystal System Corp., Japan) 
with four 300 W halogen lamps as heat sources. The growing conditions were as follows: 
The growth rate was 1 mm/h, the feeding and seeding rods were rotated at about 15 rpm 
in opposite directions to ensure the liquid’s homogeneity, and an oxygen and argon mixture 
at 3 bar pressure was applied during growth. The as-grown crystals were post-annealed at 
850°C in order to release the internal stress and to adjust the oxygen content. 

Experimental Setup 

Ultrafast reflectivity measurements were performed using a Ti;Saph amplihed laser system 
providing sub-50 fs pulses at 1.55 eV with a repetition rate of 6 kHz. The reflectivity change 
of the sample was measured in a pump-probe geometry. Of the ampliher output, 0.5 pJ 
per pulse were used to excite the sample along the c-axis; the temperature was controlled 
between 10 K and 300 K in a closed cycle cryostat. In the probe branch, the beam was 
delayed in time using a mechanical delay stage and focused into a calcium fluoride (CaF 2 ) 
crystal, generating a broadband continuum between about 1.7 and 2.9 eV. This broadband 
pulse was focused onto the sample using parabolic mirrors and its reflection collected via 
optical hber into a 4/-spectrometer. The latter used a CMOS array synchronized to the 
incoming laser pulses, enabling us to record the full probe spectrum for every single laser 
shot. Every second pump pulse was blocked by a mechanical chopper, and for each time 
delay the extracted signal is the relative change in reflectivity, AR/R, given by the difference 
of the probe spectra with and without the pump pulse. For both pump and probe beam, 
the polarization was controlled using a retardation plate. 
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The setup has an intrinsic noise level of about 0.1% RMS. The signal-to-noise ratio is 
improved statistically by repeating each measurement many times as well as binning the 
signal of several detector pixels, thus averaging in energy. In this way, relative reflectivity 
changes below 10“^ can be observed. All broadband data presented in this paper have been 
numerically corrected for the group velocity dispersion experienced by the probe pulse. 


Singular Value Decomposition Algorithm 


We present a general algorithm to separate differential reflectivity measured as a function 
of time delay and energy (^) {t, E) into its different physical components. We start by 
invoking the large separation of energy (or time) scales of the phenomena observed in real 
time {hw < 0.1 eV) and the probed window {hw > 1.5 eV) and in the spirit of the BornOp- 
penheimer approximation we assume that the differential reflectivity at high energies can 
be taken as a parametric function of a set of “slow” variables describing the out-of- 
equilibrium state produced by the pump. Here, t represents the time delay from the arrival 
of the pump pulse. Examples of slow variables are ionic displacements and “slow” charge 
and magnetic fluctuations, either coherent or incoherent. “Fast” fluctuations produced by 
the pump are either relaxed after the pump passage (~ 50 fs) or not resolved. We can thus 
expand the differential reflectivity as a function of the high probe energy E and time delay 
t as 



(E) m 


( 6 ) 


Typically, the sum can be restricted to a small number Np of processes which contribute 
signihcantly. In particular, the sum can be partitioned into an oscillatory part (t, E) 

and a relaxational part (f, E). Eq. ([^ allows to represent the full two-dimensional 

data in terms of a few time dependencies of the excitation coordinates ^i{t) and their as¬ 
sociated energy dependencies {E) which can be related to the Raman prohle for 

excitation ^ as in Eq. (|^. 

In principle one can obtain the decomposition by fitting the experimental data with a 
model of the form of Eq. (|^, but in practice it is more convenient to use the following 
algorithm. It is natural to consider the data in Fig. j^a) as a rectangular matrix (^) (t, E) 
of differential reflectivity values with the discrete version of the variables E and t playing the 
role of indices. Singular Value Decomposition (SVD) utilizes the fact that any rectangular 
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matrix (^) (t, E) G M(m, n) (corresponding to data at m time points and n energy values) 
can be uniquely decomposed into a sum of tensorial vector products [3lj of the form 


/A7?\ 

( — ((, B) = 0 v,{E) = J2 «i(() ® (7) 

^ i=l i=l 

The Aj G M are the sorted singular values, Ai > A 2 > A 3 > ..., the ® denotes the outer 
product, N is called the rank of the matrix (^) and Ui{t) G M(m) and Vi{E) G M(n) will 
be called the canonical time and energy traces, respectively. In the second form we have 
absorbed the singular value in the dehnition of the canonical time trace By construction, 
the canonical traces are orthogonal, n* • Vj = 6ij and Ui ■ Uj = Restricting the sum over 

i in Eq. ([^, one hnds that only a small number of terms are necessary to reproduce the 
physically relevant signal in (^) while the rest contribute to the background noise. We can 
therefore limit the sum to i = 1,.., W N. Note that the Ui{t), Vi{E) do not correspond 
directly to the physical quantities in Eq. 0 - This is because the canonical traces are by 
construction orthogonal while the “physical traces” j {E), ^i(t) do not need to be 

so. The latter typically consist of incoherent charge relaxations and/or damped phonon 
oscillations. We thus decompose each of the W canonical time traces Ui{t) into Np physical 
traces by htting with a sum of relaxations and damped oscillations representing the physical 
traces Ui{t) oc 

Np 

Ui{t) = 

i=i 

Inserting this into Eq. ([^ we obtain the decomposition of our data in terms of physical 
traces, 

Np 

{t,E)=j^Uj{t)e>Vj(E), (8) 

i=i 

with Vj{E) = 'Yl!i=i'^i{E)aij. We normalize ~ ^ which implies unity norm for the 

physical energy traces, Vj- ■ = 1 with Vj{E) oc {E). 

Eq. (|^ can be separated into an oscillatory part and a relaxation part, according to the 
character of the model functions Uj{t), 
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(a) (b) 




Time Delay [fs] 


Time Delay [fs] 


FIG. 4. Canonical time traces (a) ui{t) and (b) U 2 {t) obtained from the singular value decompo¬ 
sition. The traces have been fitted with a sum of two exponentials and a damped oscillation. 


For the present data we find that the oscillation of interest is present in the first two 
canonical traces Ui{t). Thus, we choose Nc = 2. We find that Ui and U 2 can be fitted 
accurately with the sum of two exponential decays (ui) or an exponential decay and a 
step funciton (M 2 ), plus one oscillation, as shown in Fig. Leaving the frequency of the 
oscillation to vary independently in the two fits, we find almost the same frequency a; = 21 
meV (169 cm“^) and a; = 20 meV (161 cm“^) respectively confirming the presence of the 
same physical trace in the two leading canonical traces. Next, we repeat the fit constraining 
the frequency and the relaxation time to be the same in the oscillatory part so that we can 
obtain the “weight” of the oscillatory physical trace in each of the two canonical traces. 
To get an accurate fit we need to introduce introduce a slight delay (18 fs) of the trace 
Ml with respect to the trace M 2 . We attribute this is to a numerical error caused by small 
uncertainties in the correction for group velocity dispersion in the original data. Thus, we 
neglect this small difference in the analysis. 

The physical time traces Ui{t) are shown in Fig. |^a), where the relaxation is contained 
in Ui to f/ 4 , whereas the oscillation is given by f/ 5 . Figure [^b) shows the corresponding 
energy vectors Vi{E). The energy dependence of the oscillation displayed in Fig. ga) is 
determined by the energy vector V 5 associated with the oscillation. The reconstruction of the 
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Time Delay [fs] 


(b) 



FIG. 5. (a),(b) Physical traces Ui{t), Vi{E) extracted from the singular value decomposition. The 
relaxation is included in Ui to 1 / 4 , while the oscillation is represented by U^. The energy vectors 
Vi have been smoothed using a sliding average. 



Time Delay [fs] 


Time Delay [fs] 
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FIG. 6. Reconstruction of the data using the SVD results up to second rank: (a) Reconstructed 
data, compare to Fig. j^a). (b), (c) Relaxation and oscillatory contributions according to Eq. (§. 


data from the SVD for = 2 and its separation into relaxation and oscillatory components 
is shown in Fig. As mentioned, the main features of the data are reconstructed even 
when limiting the SVD to second rank. Compared to the simpler and more commonly used 
method of directly using a FFT on the traces of Fig. [^b), the SVD is a lot more effective in 
subtracting the incoherent peak, reduces the noise level, and delivers the energy dependence 
of the oscillation with ease. 
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Assignment of the coherent oscillation 


The estimates for the oscillation frequency mentioned above allow to identify the coherent 
oscillation present in our data as the Aig phonon mode at 156 cm“^ present in static Raman 
data reported by Sugai for the orthorhombic phase of LCO [35]. Previous computations by 
Mostoller et al. (Table V in [HE]) show that this mode, listed as 148 cm“^ at the T point in 
the orthorhombic phase, corresponds to an Ag mode of 135 cm~^ (4 THz) at the X point 
in the tetragonal phase of LCO. The eigenvector of the latter was calculated by Falter and 
Schnetgoke (Fig. 7 of [25]): It is a breathing mode of the La atoms along the c-axis, where 
the apical O atoms move in phase with the La atoms. This mode thus differs from the 
coherently excited Aig mode in the study of optimally doped La 2 _ 3 ;Sr 3 ;Cu 04 by Mansart et 
al. [S], in which the La and apical O atoms are breathing in antiphase. We will not consider 
the minimal displacements of the in-plane O atoms indicated in [25] . 


Estimate of the electron-phonon matrix element from a point charge approximation 


The ionic displacement induced by the coherent phonon along the c-axis will affect the 
Madelung energies Scu and Sq of the copper and oxygen atoms in the ab-plane, changing the 
charge-transfer energy A = ecu — £o = Eg [57], as well as the density of states D{u) in the 
polaron band W. The electron-phonon matrix element reflecting the change in Madelung 
energies due to the displacement of ions along the c direction can be easily estimated in point 
charge approximation. To this end, we consider the effect of the displacement 5z of the La 
atoms closest to the CUO 4 plaque on the in-plane CT energy A. The atoms in question are 
the ones sitting above and below the plaque on the positions indicated in Fig. |^d). The 
shortest distances to the in-plane atoms are indicated by the dashed lines. 

The energy A is given by difference in Madelung energies of the Cu and O atoms, A = 
^Cu — = Eg [57]. Therefore, A is given as a function of the shortest La distance z from 

the Cu-0 plane by 


A(^) 


/ +4 ^ -8 \ 

ATieoER ydia-Cu dia-O J 


Ze2 

At^EoEr 


4 



- 1 - z‘^ 


8 



+ z'^ 



(10) 
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with the La charge Z = 3 and the unit cell length a = 3.8 A, which is equivalent to twice 
the in-plane Cu-0 distance. Note that there are two La atoms closest to each CUO 4 plaque 
which are nearest neighbors to two Cu atoms and four O atoms each and that the Cu and O 


atoms carry a single charge, yielding the factors four and eight in Eq. (10). For the relative 
dielectric constant Er we use the value ei{EcT) = 5.3 at the CT energy (see Fig. [Ja)) to 
account for screening effects. Using the equilibrium distance zo = 1.79 A, we find an the 
el-ph matrix element of = ^ = 4.8 eV A“L 
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